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ABSTRACT 



Thermal images of a ship target were recorded with an AGA 
Thermovision 780 thermal imaging system on 6-13 May 1987. 
These images were used to produce temperature distributions 
of the ship. The temperatures predicted with the AGA 
computer-processing algorithm were compared to actual 
temperatures measured with thermistors at eight locations on 
the ship's superstructure. An empirical modification based 
on the identity of the ambient and local atmospheric 
conditions was developed for the atmospheric compensation 
algorithm. The predicted temperatures were found to agree 
with the actual temperatures within 1.5 °C in 77 % of the 
measurements. The Contrast Transfer Function (CTF), 
Modulation Transfer Function (MTF) , and Minimum Resolvable 
Temperature Difference ( MRTD ) were determined for the 
Thermovision using flat black painted aluminum bar targets. 
The resultant curves showed the expected form with some 
experimental scatter at higher spatial frequencies. 
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I . INTRODUCTION 



Thermal imaging is the process of collecting thermal 
radiation emitted by an object and producing a visible image 
of that object. Although the concept is simple, 
implementation of a thermal imaging device quickly becomes a 
complicated venture. First, the designer is confronted with 
a relatively broad portion of the electromagnetic spectrum, 
the infrared region; therefore, he must consider target 
radiation characteristics in the design process. Second, the 
target must be recognized against a background of several 
natural radiation sources. These sources include celestial 
emissions, the sun, the sky, and terrestrial materials. 
Finally, the atmosphere attenuates target radiation by 
scattering and absorption mechanisms. Most current thermal 
imagers utilize quantum detectors to collect infrared 
radiation from a scene and convert this information to analog 
electrical signals which are processed and displayed on a 
video monitor. The ability of a detector to recognize a 
target depends on image contrast, i.e. the difference between 
target emittance and the background emittance. Because 
thermal radiation is a function of an emitter's temperature, 
the effective temperature difference between target and 
background provides a good measure of image contrast provided 
the emissivity differences are not great. 
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The apparent temperature perceived by an imaging device 
can vary considerably from the actual temperature. This 
irregularity is primarily due to background and atmospheric 
effects. If a designer is able to predict the impact of 
these effects accurately, then the imaging device will be 
better able to match perceived target temperatures with 
actual target temperatures. 

This thesis explored the differences between the apparent 
temperatures sensed by a thermal imaging system and the 
actual temperatures of a target. The imaging system used was 
an AGA Thermovision 780, manufactured by the AGEMA 
Corporation. The target was the R/V Point Sur, an 
oceanographic ship operated by the National Science 
Foundation for the Naval Postgraduate School and Moss Landing 
Marine Laboratory. The background was Monterey Bay in two 
vicinities: the Moss Landing Marine Laboratory and the 
Hopkins Marine Station. Data were taken during May 1987. 
Also, the Thermovision's Modulation Transfer Function ( MTF ) 
and Minimum Resolvable Temperature Difference (MRTD) curves 
were determined using square aluminum plates painted with 
bars of flat black paint. The use of aluminum permitted 
adequate heat transfer across a plate, so that the 
temperature was constant at all portions of the plate. The 
emissivities of the paint and aluminum were substantially 
different which caused a radiation contrast to exist between 
the bars of paint and bars of exposed aluminum. 
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Chapter II details the theory behind thermal radiation, 
atmospheric propagation of such radiation, and thermal 
imaging system performance measures. Chapter III describes 
the technical characteristics and operational capabilities of 
the Thermovision 780, as well as the temperature measurement 
techniques of the system. Chapter IV contains the results of 
the experiment, while the analysis of these results is 
presented in Chapter V. The conclusions and recommendations 
of this thesis follow in Chapter VI. 
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II. PERTINENT THEORY 



A. THERMAL RADIATION LAWS 

Thermal radiation is comprised of photons that are 
emitted by molecules undergoing vibrational and rotational 
transitions. These molecules may be the original sources of 
the photons or they may be reradiating photons emitted from 
other sources. In macroscopic terms, objects can be self- 
emitting thermal radiation sources or reflectors of such 
radiation. Most objects in our environment exhibit both 
mechanisms. [Ref. l:p. 18] 

Several physical laws have been developed that describe 
the nature of thermal radiation. Foremost is Planck's 
radiation law which gives the spectral distribution of 
radiant emittance, W^. This law is: 

W^ [ X , T ] = (2nhc/X 5 ) 1 (Watt/m 2 - pm) (2-1) 

[exp(ch/XkT) - 1 ] 

where k is Boltzmann's constant, c is the speed of light in a 
vacuum, and h is Planck's constant. Figure 2.1 [Ref. 1: 
p. 23] depicts spectral emittance for three object 
temperatures. When this expression is differentiated and set 
equal to zero, Wien's displacement law is obtained. This is: 

X m T = 2897.8 {pm-K} (2-2) 

where X m is the most abundant wavelength of radiation for a 
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particular temperature. This thesis deals with the 8-14 
micrometer band of radiation which corresponds to a 
temperature range of 207.0-362.2 K, temperatures that clearly 
exist in our environment. This is a good indication of the 
background problem that confronts thermal imaging systems. 

Integrating the expression for Planck's law across the 
entire spectrum of wavelengths yields the Stefan-Boltzmann 
law, given by: 

W(T) = ( 2Tt 5 k 4 ) / ( 15c 2 h 3 ) T 4 = oT 4 {Watt/m 2 } (2-3) 

where W(T) is radiant emittance and o is the Stefan-Boltzmann 
constant. Radiation sources are commonly described in terms 
of radiance, L(T), which is radiant emittance through a unit 
solid angle. Many sources effectively radiate into a 
hemisphere as perfectly rough planes. Radiance for such a 
Lambertian surface is: 

L(T) = W ( T ) / it {watt/m 2 -steradian} (2-4) 

Most objects in our environment do not emit as ideal 
blackbodies, but emit only a fraction of the blackbody 
radiant power. This fraction is known as the "emissivity" 
and can range in value from zero (a perfect reflector) to one 
(a black body). Thus, the Stefan-Boltzmann law becomes: 

W ( T ) = EOT 4 (2-5) 
where e is the emissivity of an object. For an object in 
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thermal equilibrium with its surroundings, the object's 
absorbed power must equal the emitted power. This leads to 
Kirchoff's law which is: 

e(M = a(\) (2-6) 

where a(\) is the object's spectral absorptivity. 

Thermal radiation impinging on an object can undergo one 
of three mechanisms. First, the radiation can be absorbed 
and subsequently reradiated as described by Kirchoff's law. 
Second, the radiation can be reflected off the object. 
Third, the radiation can be transmitted through the object. 
Hence, the boundary conditions at an object's surface are 
described by the total power law: 

a(\) + p(\) + T(\) = 1 ( 2-7 ) 

where p(M is the spectral reflectivity and T(\) is the 
spectral transmissivity. 

B. NATURAL RADIATION SOURCES 

Celestial background radiation includes infrared emission 
from solid material within our galaxy, as well as from 
extragalactic sources. Irradiance levels from celestial 
sources have been found to be around 10 -16 W/cm J -ym for the 3 
to 26 micrometer region [Ref. 2:p. 3-28]. Such levels are 

several orders of magnitude lower than those of other 
sources; therefore, celestial sources will not be considered 
in this thesis. 
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The sun can be treated approximately as a blackbody at 
5900 K, when the earth's atmospheric effects are ignored. 
Figure 2.2 [Ref. 2:p. 3-34] depicts the spectral irradiance 
of the sun outside and within the earth's atmosphere. At 
5900 K the most abundant wavelengths radiating from the sun 
are in the visible region; however, a small percentage of the 
sun's spectral emittance is in the infrared region i.e., 
0.0986 percent in the 8-14 micrometer band, which equates to 
133.5 W/m 2 [Ref. 2:p.3-36]. Contributions of the sun to 
background radiation must be considered during daylight hours 
and particularly in background scenes that approach the sun's 
direction. 

The sky contributes to the background radiation through 
the scattering of solar radiation and emission from 
atmospheric constituents [Ref. 2:p. 3-71]. Figure 2.3 [Ref. 
2 : p . 3-71] depicts these contributions as well as solar 
radiation and radiation from a cloud. The presence of clouds 
affects both the scattering of near-infrared solar radiation 
and the sky's thermal radiation. Forward scattering in 
clouds is prevalent in the near-infrared region; however, 
multiple scattering dominates the forward scattering effects 
in a heavily overcast sky. The spectral emittance due to the 
sky emissions approaches that of a black body at atmospheric 
temperature in bands of high atmospheric absorption. 

Optically thick clouds are also well represented by a 
black body curve based on a cloud's average temperature; 
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however, the "wings" of this curve correspond to regions of 
high atmospheric absorption. Accordingly, the black body 
curve in these "wings" is based on the atmospheric 
temperature rather than the cloud's temperature. Figure 2.4 
[Ref. 2:p. 3-73] illustrates this effect of a cumulus cloud 
radiating as a blackbody at one temperature and the sky 
radiating as a black body at different temperature. 

While the radiation of numerous terrestrial materials 
could be examined, only that of marine backgrounds will be 
considered in this thesis. The marine background is affected 
by four factors: the optical properties of water; the surface 
geometry and wave slope distribution; surface temperature 
distribution; and bottom material properties [Ref. 2:p. 3- 
105] . 

Sea water has very low transmittance and reflectance 
values in the infrared region. The exception is that high 
reflectances do exist at large angles of incidence. When the 
sea surface is roughened by wind, the reflectance is 
significantly reduced at angles of incidence near the 
horizon. Figure 2.5 [Ref. 2:p. 3-106] depicts reflectivity 
and emissivity of smooth water versus angle of incidence. 

As with any thermal emitter, the sea surface temperature 
determines the radiant emittance of the sea. The temperature 
gradient is quite pronounced in the upper 1.0 mm of the 
surface due to the cooling effects of evaporation. Figure 
2.6 [Ref. 2:p. 3-109] depicts this gradient. The gradient is 
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much less at lower depths but can be influenced by the 
convective activity of the water. Surface temperatures can 
be altered by surface contamination which restricts the flow 
of heat from lower depths; thus the apparent surface 
temperature will be lower than the bulk water temperature 
[Ref. 2 : p . 3-109]. 

Because sea water has such low transmittance values for 
the infrared region, the properties of bottom materials have 
negligible effects on the background radiation. 

C. ATMOSPHERIC PROPAGATION OF THERMAL RADIATION 

The transmittance of the earth's atmosphere is less than 
one, since the atmosphere does not behave as a perfect 
dielectric. This transmittance is given by the Lambert-Beer 
law : 



where T A (X) is the extinction coefficient and R is the 
distance of propagation. The average transmittance for a 
particular bandwidth X^ to X 2 i- s: 



Extinction of thermal radiation is caused by absorption and 
scattering processes. Thus, the extinction coefficient is 
given by the expression: 



T A (X) = exp[-p(X)R] 



( 2 - 8 ) 




(2-9) 



M(X) = k (X) + a ( X ) 



( 2 - 10 ) 
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where k(X) is the absorption coefficient and a(X) is the 
scattering coefficient. Both coefficients have components 
due to air molecules, as well as aerosol particles. 
Therefore : 

k(X) = kpp ( X ) + kg(X) (2—11) 



and 



a ( X ) = a m (X) + a a ( X) (2-12) 

where the subscripts m and a denote molecule and aerosol 
respectively. The values of these coefficients have been 
empirically determined and depend on the density and 
composition of both the molecules and aerosols. 

Molecular absorption is due primarily to water, carbon 
dioxide, ozone, nitrous oxide, carbon monoxide, and methane 
[Ref. 2:pp. 5-101 - 5-105]. These absorbers limit thermal 

radiation in the atmosphere to two windows: the 3. 5-5.0 and 

8-14 micrometer bands. 

Aerosols have been categorized in four standard types. 
Maritime aerosols are made up primarily of salt particles 
which act as condensation centers for water. The 

concentration of these particles is largely dependent on wind 
speed, while the size distribution is influenced by relative 
humidity. Continental aerosols consist of silicon, iron, 
sulphates, and organic material. One third of these aerosols 
can act as condensation centers. Urban aerosols are composed 
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of combustion and industrial products. Stratospheric 
aerosols contain sulphate with an occasional addition of 
volcanic dust. [Ref. 3: p. 17-14] 

Scattering by molecules in the 8-14 micrometer band is 
relatively insignificant. Wavelengths in this region are 
much larger than particle size, so that Rayleigh scattering 
would be valid. However, the X -4 dependence of Rayleigh 
scattering reduces the effect of molecular scattering. Mie 
scattering theory pertains to particles of all sizes; 
therefore, it is valid for aerosols. Because Mie scattering 
is somewhat independent of wavelength in the large particle 
limit, aerosol scattering is quite significant in the 8-14 
micrometer band. Mie scattering is most effective around 
r/\ - l, while particle size distributions extend out to tens 
of micrometers. 

Due to the complex nature of atmospheric propagation of 
thermal radiation, several computer models have been 
developed to predict atmospheric transmittance. LOWTRAN 6 is 
the model which will be used in this thesis. LOWTRAN is a 
FORTRAN based code that calculates atmospheric transmittance 
and radiance for a specified range of wavelengths along a 
designated path length. Seven atmospheric models are 
available including the U.S. Standard Atmosphere and a model 
based on user specified meteorological data. 
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D. PERFORMANCE MEASURES OF THERMAL IMAGING SYSTEMS 

A fundamental parameter restricting performance of any 
imaging system is contrast. In the case of thermal imagers, 
radiation contrast is used to assess how well a target can be 
seen against its background. The expression for radiation 
contrast, C , is: 

C = (Wrp - Wg ) / (Wrp + Wg) (2-13) 

where w T is the target radiant emittance and W B is background 
radiant emittance. Figure 2.7 [Ref. l:p. 29] shows radiation 
contrast curves for four background temperatures. 

The radiation contrast of bar targets can be used to 
construct a Contrast Transfer Function (CTF) . As a bar 
target represents a square wave input to an imaging system, 
the CTF describes the system's square wave amplitude response 
at the spatial frequency of that target [Ref. 4:p. 114]. The 
CTF of a system is determined by plotting the radiation 
contrast for a range of spatial frequencies. The resultant 
function has an initial value of 1.0 (100 percent contrast) 
at zero spatial frequency and drops to a final value of 0.0 
(no contrast) at the system's cutoff frequency, i.e. the 
maximum spatial frequency above which a system can no longer 
resolve a test target. 

A more useful measure of a system's resolution capability 
is the Modulation Transfer Function (MTF) . "The MTF is the 
sine-wave spatial frequency response" [Ref. 4:p. 114]. The 
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MTF can be determined from the CTF using the following 
relationship [Ref. 4:p.ll7]: 

M ( f ) - ( tt/4 ) [C( f ) + C( 3f ) / 3 - C(5f)/5 + C(7f)/7] (2-14), 

where M(f) is the value of the MTF at a single spatial 
frequency, f. C(f) is the value of the CTF at that 
frequency. Figure 2.8 [Ref. l:p. 19] contains an MTF curve 
for an example imaging system. 

While the MTF is a valuable measure, it describes only a 
small part of an imaging system's ability to perform intended 
tasks. Ideally, a performance measure must be based on 
fundamental system parameters and "must relate to the 
performance of the system as it is intended . " The most 

widely used performance measure for infrared detection 
systems is Noise Equivalent Temperature Difference (NETD). 
NETD is a measure of a system's ability to detect small 
signals in noise. NETD is the target-to-background 

temperature difference that corresponds to a system's peak 
signal to rms noise ratio equal to one. [Ref. l:p. 166] 

NETD is derived from an expression for spectral irradiant 
power received by the detector. Such an expression would be: 

P\(\,T) = [W x (\,T)/tt] 0A t T o (Watt/ym) (2-15) 

where W^(\,T) is the target's spectral emmittance, 0 is the 
solid angle subtended by the system's optical aperture, A T is 
the target area from which radiation is received by the 
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system's optical aperature, and T 0 is the system optical 
transmission. Since the solid angle is A 0 /R 2 and the target 
area is apR 2 , equation 2-15 can be written as: 

P\ = (Wx/tt) A 0 a(3T 0 (2-16) 

where A 0 is the area of the optical aperature and a and p are 
the system's horizontal and vertical subtense angles. 

The target's differential change in irradiant power with 
respect to temperature is of primary interest. By 
differentiating equation 2-16 with respect to temperature, 
such an expression is found: 

8P X /8T = [ (a(3A 0 T 0 )/TT] 8W X /8T (2-17) 

The system's differential change in signal voltage is 
found by multiplying equation 2-17 by the system's 
responsivity , R(X), a parameter that gives the ratio of 
signal voltage output to incident power: 

8V S /3T = R ( X ) [ ( aPA 0 T 0 ) /n ] 8W X /8T (2-18) 

Responsivity is given by: 

R(M = [V n D*(M] / ( abB n ) 1/2 (2-19) 

where V n is the detector rms noise voltage produced by B n , 
the noise bandwidth of a test reference filter, a and b are 
detector dimensions, and D*(\) is the system's specific 
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detectivity. Inserting equation 2-19 into 2-18 yields: 

8V S /8T = D*(A) (dPA 0 T 0 ) / [ tt ( abB n ) 1 1 2 ] 8W X /8T (2-20) 

Assuming small target-to-background temperatures, the 
temperature differential of target emittance can be 
approximated by the derivative of Planck's law at the 
background temperature T B : 

8W X /8T « (C 2 /VT B 2 ) W X (T B ) (2-21) 

where c 2 = 1.4388 X 10 4 ym-K. Specific detectivity can also 
be written as: 



D*(\) = AD*(A p )/A p 
= 0 



for X S Aj 
for A 2 X T 



( 2-22 



where Ap is the wavelength corresponding to the peak value of 
system detectivity. Inserting equations 2-21 and 2-22 into 
equation 2-20 and integrating across the effective waveband 
A]_ to Ap yields: 

8V S /8T = o^AoToYn D* ( A p ) c 2 f X P W X (T B ) dA (2-23) 

tt ( abB n ) 1/2 Ap T b 2 J X1 



Assuming a small signal approximation and rearranging 
equation 2-23 produces an expression for the signal-to-noise 
ratio : 



SNR = 



^ v s/ v n 



= AT 



0PA 0 T 0 



D*(X p )c ? f Xp W X (T B ) dA (2-24) 



TT(^bB-)l/ 2 — A^ J X1 
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Based on the definition of NETD, the SNR is set equal to 
one and AT becomes NETD. Thus, NETD in final form for a 
scanning imager is: 



Although NETD is a satisfactory measure of a system's 
target detection capability, it is of limited utility for 
thermal imagers because this measure does not account for 
image quality. A more appropriate performance measure for 
these devices is Minimum Resolvable Temperature Difference 
(MRTD) which incorporates the system MTF and provides a 
measure of how well an imaging system can resolve a target. 
MRTD is defined as the blackbody target-to-background 
temperature difference in a standard test pattern at which an 
observer viewing the system display can resolve a target [Ref 
1 : pp . 190-192]. 

An analytical expression for MRTD is based on a system's 
signal-to-noise ratio in the image of one bar. This ratio 
is : 



where M(f) is the system modulation of the image (i.e. MTF), 
the AT/NETD factor is the electronic SNR measured at a NETD 
reference filter with a target-to-background temperature 




SNRj_ = 4/tt M ( f ) [AT / (NETD] (p 1/2 ) _1 



(2-26) 



difference AT, and p 1/2 is the ratio of actual system 



bandwidth to reference bandwidth. 
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As the observor views the target, this perceived SNR is 
modified by four factors. First, the eye operates on a mean 
signal rather than peak signal; therefore, the SNR is reduced 
by a factor of 2/tt. Second, the SNR is improved by a 
temporal integration factor of (TgF)- 1 -/ 2 where T e is the 
effective eye integration time and F is the system frame 
rate. Third, the SNR is improved by the eye's spatial 
integration of the bar height which is seven times greater 
than the bar width. This factor is (7/2 f (3 ) 1 ^ 2 . Finally, 
the bandwidth ratio is improved by the eye's matched filter 
action. [Ref. l:p. 186] 

A good approximate value of the perceived SNR has been 
found to be 4.5 for a 90% probability of detection [Ref. l:p. 
188]. Using this value and the observer's perception 

factors, equation 2-26 is rearranged so that: 

MRTD(f) = AT = 3 NETD p m 1/2 (fp) 1 / 2 (2-27) 

M ( f ) (T e F) 1/2 

where Pm 1 ^ 2 is the improved bandwidth ratio due to the eye's 
matched filter action. Because the displayed noise in many 
systems is white [Ref. l:p. 189], p^^ 2 can be given by the 
simple expression: 

p m 1/2 = ( 4af/TT) 1/2 (2-28) 

Inserting equation 2-28 into 2-27, the final expression 
for MRTD becomes: 



23 



MRTD(f) = 6 NETD f(a&) 1/2 
M(f) (nT e F) 1/2 



( 2-29 ) 



{K} 

Figure 2.8 contains an MRTD curve for an example imaging 
system. 

The NETD listed for the AGA Thermovision 780 is 0.12 °C 
at 22 °C [Ref. 5]. The MRTD is not provided by the 
manufacturer but was determined during the experimental 
portion of this thesis. These results are found in Chapter 
IV. The MRTD deduced is then applied to the discussion of 
the measured temperature distributions in Chapter V. 
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III. DESCRIPTION OF AGA THERMOVISION 780 

A. GENERAL 

The AGA Thermovision 780 thermal imaging system 

incorporates a dual scanner which senses thermal radiation in 
two spectral bands and produces electronic video signals. 
These signals are amplified and used to display images on a 
black and white monitor. The scanner and monitor are shown 
in Figure 3.1 [Ref. 6]. The system is augmented with a 
microcomputer that displays digitized pictures of the 

system's video on a color monitor and records these pictures 
onto a floppy disk. The computer uses a software program 
that color codes each picture according to the scene's 
temperature distribution. A schematic of the total 

Thermovision system is displayed in Figure 3.2 [Ref. 7]. 

B . DUAL SCANNER 

The dual scanner is actually comprised of two scanning 
systems. The shortwave system, 3 to 5 . 6 micrometer, has a 
single Indium Antimonide photovoltaic detector and silicon 
lens. The longwave system, 8 to 14 micrometer, has a single 
Mercury Cadmium Telluride detector and germanium lens. Each 
system has a 7° X 7° field of view with a standard lens. A 
3.5° X 3.5° lens is also available for the long wave system. 
Each system uses vertical and horizontal scanning prisms to 
produce a 4 to 1 interlaced raster scan. Each of the four 
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fields has 100 scanning lines, although only 70 lines are 
used to produce images. The interlaced frame in this 
configuration consists of 280 scanning lines. The detectors 
are cooled to 77 K with Dewar flasks containing liquid 
nitrogen . 

C. BLACK AND WHITE MONITOR CHASSIS 

The dual scanner is connected to two black and white 
monitors, so that scenes in the two spectral bands can be 
viewed simultaneously. The monitor chassis contains controls 
for adjusting the brightness and contrast of the display. 
The chassis also contains controls for adjusting the thermal 
level and thermal range of the system. These two adjustments 
are measured in Thermal Units which are arbitrary units of 
measure proportional to the intensity of the system's 
received thermal radiation. The thermal level control 
adjusts the DC level of the AC video signal, while the 
thermal range control limits the dynamic range of this 
signal . 

These controls can be used in conjunction with a 
calibration plot to assess manually the temperature 
distribution of a scene. The relationship between thermal 
level and thermal range settings versus temperature is shown 
in the calibration curve enclosed as Figure 3.3 [Ref. 6]. A 
change in thermal level corresponds to a nonlinear change in 
scene temperatures the Thermovision can sense. An adjustment 
of the thermal range determines the measurable temperature 
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range which is approximately centered about the "median" 
temperature established by the thermal level setting. 

D. IF 800 MICROCOMPUTER 

A BMC IF 200 microcomputer processes the infrared data 
for the AGA Thermovision 780 and provides an automated means 
of assessing temperature distributions. On command, the 
computer's DISCO 3.0 program (AGEMA Corporation proprietary 
software) creates an eight color picture based on thermal 
values provided by the system and parameters inserted by the 
user. The color scheme depicts the temperature distribution 
of a scene with each color representing a particular 
temperature range. The program provides temperature data to 
the user in several formats. The feature utilized in this 
thesis is the ability to assess the temperature associated 
with each pixel in the display screen comprised of an array 
of 128 X 128 pixels. 

The computer in this temperature evaluation mode displays 
a crosshair which is moved about the screen by means of four 
cursor control keys. The user positions the crosshair over a 
particular pixel and the computer displays the pixel's 
temperature in the lower right corner of the screen. 

E. TEMPERATURE MEASUREMENT TECHNIQUE 

The thermal measurement technique utilized by the 
Thermovision 780 is based on the relation: 

p i - T a £ o p o + T a (l-e 0 )P s + (l-T a )P a {Watt} (3-1) 



27 



where Pj_ is the total radiant power received by the system, 
P Q is the radiant power from the object as a blackbody, P s is 
the radiant power from the object's surroundings as a 
blackbody, P a is the radiant power from the atmosphere as a 
blackbody, T a is the atmospheric transmittance, and e 0 is the 
object's emissivity. The emissivity of the surroundings and 
atmosphere are assumed to be one. 

The first term on the right hand side of equation 3-1 is 
the received radiant power emitted by the object. The second 
term is the received radiant power emitted by the object's 
surroundings and reflected by the object. The third term is 
received radiant power emitted by the atmosphere. 

Because the system's thermal value measure is 
proportional to received radiant power, equation 3-1 can be 
written as: 

li = T a s 0 I 0 T a ( 1-e o ) J s (1 — T a )la (Thermal Units} (3-2) 

where the I terms represent the thermal values of 
corresponding radiation sources. The received thermal value 
consists of two terms such that: 

Ij_ = L + i ( 3-3 ) 

where L is the thermal level setting on the monitor chassis 
and i is a fractional portion of the thermal range. 

Substituting equation 3-3 into 3-2 and manipulating terms 
yields the following expression for the object's thermal 
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value : 



I Q = (L+i)/T a e 0 - (l/e 0 -l)I s - l/e Q ( 1/T a -1 ) I a (3-4) 

The thermal values in the above expression are dependent on 
corresponding temperatures. The thermal value-temperature 
relationship is calibrated under laboratory conditions using 
the equation: 

I = A / [C exp ( B/T ) -1 ] (3-5) 

where A,B, and C are calibration constants. The resultant 
calibration curves are similar to that shown in Figure 3.3. 
Equation 3-5 is used by the system to convert atmospheric and 
ambient temperatures to corresponding thermal values and to 
convert I to an object temperature. 

The atmospheric transmission factor is approximated by 
the system using: 

T a = exp[-a(/d-l ) ] ( 3-6 ) 

where a is an atmospheric attenuation constant and d is the 
distance to object. The value of a specified by the 
manufacturer is 0.008 for the long waveband. A more accurate 
transmission factor can be computed using the LOWTRAN 
propagation code which is based on a form of Beer's law as 
given in equation 2-8: 

T a (\) = exp[ -y ( X ) R ] (3-7) 
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where y(\) is the atmospheric extinction factor that 
represents scattering and absorption effects as discussed in 
Chapter II. Values based on this equation will be used in 
this thesis. 
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IV . ACQUIRED DATA 



A. RECORDING OF SHIP IMAGES 

The AGA Thermovision 780 was used to record thermal 
images in the 8-14 micrometer band of the R/v Point Sur at 
the Moss Landing location on 6 May and 13 May 1987 and at the 
Hopkins Marine Station on 9 and 11 May 1987. The map in 
Figure 4.1 depicts these locations on Monterey Bay and the 
approximate direction in which the Thermovision was aimed. 
Images were recorded with the ship presenting various 

aspects; however, only images of the ship at a 90° aspect 
were evaluated in this thesis. The Point Sur was 

instrumented with eight thermistors to measure the 

temperatures of various parts of the ship. These 

temperatures would be used as a basis of comparison for 
temperatures sensed by the Thermovision. Weather balloons 
were launched prior to each recording period. The radiosonde 
data recorded by the balloons were used as input to the 
LOWTRAN 6 program to determine atmospheric transmittance 
values. These data are contained in Table 1. 

Images were recorded with the Thermovision's thermal 
range setting on two, five, and 10 Thermal Units; however, 
only those images recorded at the setting of five were used 
to evaluate temperature distributions of the ship. A thermal 
range of two Thermal Units did not correspond to a large 
enough temperature range in order to measure the entire ship. 
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A thermal range of 10 Thermal Units produced images with 
poorer temperature resolution than those images recorded at a 
setting of five Thermal Units. 

B. EVALUATION OF SHIP IMAGES 

Before the recorded images could be evaluated for 

temperature distributions, the Thermovision had to be 
initialized with the following values: the ambient air 

temperature around the ship; the atmospheric temperature; the 
emissivity of the ship; the range of the ship; and the 
atmospheric transmittance. The values of ambient temperature 
and the atmospheric temperature were assumed to be equal and 
were based on the air temperatures reported by ship 

personnel. These temperatures are contained in Table 2. The 
emissivity of the ship was assumed to be 0.95 as previously 
determined and reported [Ref. 9]. The range of the ship from 
the Thermovision's location was determined by measuring the 
length of the ship's image on the Thermovision's display 
screen and inserting this value into the following 

trigonometric relation: 

R = (S Q /tan 9) (S-j/w) -1 = (wS Q /tan 9) 1/S± (m) (4-1) 

where R is the ship's range, w is the display sceen width 
(13cm), S G is the ship length (41.2m), 9 is the system field 
of view (7°), and is the image length in cm. Inserting 
values for the constants, equation 4-1 becomes: 

R = (4354 / s i ) ± 20 {m) (4-2) 
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where the uncertainty in range was based on the differential 
method of error estimation. 

As mentioned previously, the transmittance values were 
calculated with the LOWTRAN 6 code. The atmospheric model 
was based on radiosonde data collected by the weather 
balloons. The Navy maritime aerosol model was used which 
incorporated three aerosol components: a continental 
component; a "stationary" component produced by winds and 
whitecaps; and a "fresh" component formed by current 
conditions [Ref. 3 : p - 17-35]. 

An air mass character of 2 was selected for the Moss 
Landing location, while an air mass character of 4 was 
selected for the Hopkins Marine Station. The air mass 
character was a subjective rating of continental influence on 
maritime aerosols on a scale of one to 10. A rating of one 
would represent an open ocean, while ten would indicate 
strong continental influence. As can be seen in the map 
enclosed as Figure 4.1, the Moss Landing location was quite 
exposed to the open ocean, while the Hopkins Marine Station 
was somewhat protected by the Monterey Peninsula. 

A horizontal path was specified with an altitude of 5 m. 
The range was based on the predicted range for a given image. 
Calculated transmittance values are contained in Table 2. 

The Thermovision's image evaluation program also had to 
be calibrated to a reference heat source prior to evaluation 
of the Point Sur images. The calibration was conducted on 
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13 October 1987 using a black body heat source adjusted to 
temperatures ranging from 6.0 °C to 30.0 °C. The black body 
consisted of a hollow aluminum cylinder covered with fiber- 
glass insulation. The cylinder was 28 cm long and 7.5 cm in 
diameter. A 2 cm entrance hole was located in the front of 
the cylinder. The rear surface was canted, so that incident 
light normal to the entrance hole would not be reflected out 
of the cylinder. A thermocouple was attached to the rear of 
the cylinder for the purpose of measuring the black body 
temperature. The interior of the cylinder was painted flat 
black . 

The calibration was conducted in a laboratory with a room 
temperature of 18.8 8 C. In order to cool the black body to 
below room temperature, liquid nitrogen was poured into the 
cylinder until the black body temperature was 2 °C. As the 
black body warmed to room temperature, the Thermovision was 
used to record images of the black body at 1 °C increments 
starting at 6 °C. In order to warm the black body above room 
temperature , the cylinder was wrapped with a nichrome wire. 
The wire ends were connected to a variac, so that the current 
generated would heat the cylinder. Again the Thermovision 
was used to record images of the black body in 1 °C 
increments up to 30 °C. 

Once the black body images were recorded for each 
temperature, the Thermovision's image evaluation program was 
used to predict the black body temperatures. The results of 
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the calibration process are contained in Table 3. These 
results were' used to produce the calibration curve shown in 
Figure 4.2. This curve was subsequently used to convert 
temperatures computed with the Thermovision method to values 
that more closely approximated actual temperatures. 

The temperature distributions of the Point Sur were 
developed using a 6 X 15 element array. This array 
configuration was based on the Point Sur's approximate height 
to length ratio of 1/5. The vertical and horizontal 
dimensions were multiplied by factors of six and three 
respectively in order to facilitate the presentation of 
distributions on 8 1/2" X 11" notebook paper. In this 
configuration one distribution element approximately 
represented a 2.75 m X 1.35 m portion of the ship's surface 
area at 90° aspect. Images used to produce these 
distributions were recorded during the following times: 



DATE 
6 May 
6 May 
9 May 
9 May 
11 May 
11 May 
13 May 
13 May 



TIME 


SIDE OF SH 


1857 


Starboard 


1901 


Portside 


1128 


Portside 


1140-1141 


Starboard 


0944-0948 


Starboard 


1004-1008 


Portside 


1520-1521 


Starboard 


1526 


Portside 
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The temperature distributions for these periods are found 
in Figures 4.3 through 4.6. An anomoly was found during the 
initial review of these distributions. The inconsistency 
pertained to Figure 4.3 which portrayed the Point Sur's 
temperature distribution for 6 May. Some temperatures of the 
ship's upper superstructure in this figure were lower than 
the ambient air temperature. Physically, this situation 
would not be possible. The Point Sur had no known heat sinks 
in this portion of the ship; therefore, ship temperatures 
could equal or exceed the ambient air temperature, but would 
not be less than the air temperature. The conclusion drawn 
from Figure 4.3 was that the Thermovision might have been 
underestimating actual temperatures. This problem had to be 
resolved before any detailed analysis of the Point Sur's 
temperature distributions could be made. 

C. THERMISTOR DATA 

The thermistor resistance values were recorded by 
microcomputer throughout the time the Point Sur was at sea. 
The locations of the eight thermistors are indicated in 
Figure 4.7. After the cruise these values were converted to 
temperature values and averaged for the time periods 
corresponding to those periods for which the Thermovision was 
recording images. Actually, the thermistor averaging time 
was extended two minutes prior to and two minutes after the 
the image recording period. The averaged thermistor 
temperatures with standard deviations are found in Table 4. 
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These averaged temperatures were then compared with the 
temperatures- computed by the Thermovision for areas of the 
ship where the thermistors had been located. The temperature 
differences between the thermistor-measured temperatures and 
the Thermovision-computed temperatures are listed in Table 5. 
These differences were quite significant on 6 May and 13 May. 
The Thermovision underestimated the thermistor-measured 
temperatures by an average of 4.5 °C on 6 May and 4.7 °C on 
13 May. This situation reinforced the tentative conclusion 
that the Thermovision was underestimating actual 
temperatures . 



C. CTF , MTF, AND MRTD CURVES 

Because a Minimum Resolvable Temperature Difference 
(MRTD) curve was not provided by the manufacturer, one was 
produced for the purposes of this thesis. As discussed in 
Chapter II, the MRTD curve provides a good measure of a 
system's resolution capability. Since the MRTD curve is 
based on a system Modulation Transfer Function (MTF), the 
Thermovision's MTF had to be determined. A system Contrast 
Transfer Function (CTF) was more readily attainable; 
therefore, the system CTF was determined and the MTF computed 
from this curve. 

A Contrast Transfer Function (CTF) for the Thermovision 
was determined on 19 October 1987. The CTF was based on 
values of radiation contrast for targets of various spatial 
frequencies as described in Chapter II. The targets consisted 
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of 1/8 inch thick square aluminum plates with four bars of 
flat black paint applied. The emissivity of the paint was 
estimated at 0.95 and the emissivity of the aluminum was 
estimated at 0.30 [Ref. 10]. Thus, a target appeared to the 
Thermovision as a standard seven bar target as depicted in 
Figure 4.8 [Ref. l:p. 76]. Bar widths ranged in size from 5 
mm to 120 mm. The targets were placed at distances ranging 
from 1 m to 4 m. The following relation was used to generate 
28 various spatial frequencies: 

f = d / 2w (cycles/radian) (4-3) 

where f is spatial frequency, d is distance between target 
and Thermovision optics, and w is target bar width. 

With the laboratory room temperature at 18.8 °C the 
targets were heated to 25 °C. A thermocouple was bolted to 
the center of each target to measure target temperature. 
While the Thermovision was imaging a bar target, an 
oscilloscope was used to measure the radiation contrast 
between bars. The oscilloscope was connected to a video 
output port of the anologue-to-digital convertor and set to 
display target signal voltage. The intensity profile in 
Figure 4.8 depicts how a bar target would appear on the 
oscilloscope display. The heights of the signal peak and 
minimum were measured and equation 2-8 used to calculate 
radiation contrast. 

C = (W T - W B )/(W T + W B ) (2-8) 
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In this case W T was proportional to the emittance of the 
black bars of paint and W B was proportional to the emittance 
of the aluminum bars. (Note: The vertical distance on the 

oscilloscope display was proportional to signal voltage and 
thus target emittance.) Since the target temperature and 
background temperature were equal, the radiation contrast 
could be written as: 

C = ( SipOTip^ - EgOTg^ ) / ( SrpOTrp^ + EgOTrp^ ) 

= (E T - e B )/(e T + e B ) (4-4) 

The values of radiation contrast for various spatial 

frequencies are listed in Table 6. These values were 
normalized relative to 4 cycles/radian , the frequency at 
which the bar target consisted of one black bar and one 
aluminum bar. This target filled the system's entire field 
of view, hence this was the minimum spatial frequency that 
could be viewed by the Thermovision. The normalized 
radiation contrast values were then used to calculate the MTF 
values for each spatial frequency based on equation 2-14: 

M ( f ) - (TT/4)[C(f) + C ( 3f ) / 3 - C(5f)/5 + C(7f)/7] (2-14) 

CTF and MTF values are found in Table 7. Finally, the MRTD 

was calculated for each spatial frequency using equation 
2-29 : 



MRTD ( f ) 



6 



NETD f(a(3) 1/2 
M ( f ) (TTTgF) 1 ^ 



(2-29) 
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The NETD for the Thermovision was listed by the 
manufacturer as 0.12 8 C at a background temperature of 22 °C. 
However, the MRTD calculated in this thesis was based on a 
background temperature of 18.8 °C; therefore, a correction 

factor was needed. This factor was readily computed using 
the expression for NETD found in equation 2-25: 

NETD = tt ( abB n ) 1/2 X p T R 2 rf X P W X (T R ) dX^ -1 (2-25) 

a(3A 0 T 0 D*U p )c 2 LJ Xl 

The factor dependent on background temperature was: 

CF(T b ) = T b 2 -f X P W X (T B ) dX- _1 (4-5) 

J X 1 

The denominator of this factor was integrated from X^ = 8 ym 

to Xp = 10 ym with a hand-held calculator using Simpson's 
rule. The corrected NETD was: 

NETD (18.8 °C) = CF (18.8 °C) NETD ( 22 °C) = 0.985 (0.12 °C) 

CF ( 22 °C) 

- 0.12 °C 

Thus, the NETD was essentially unchanged. 

The following Thermovision parameters were inserted into 
equation 2-29: 

NETD = 0.12 °C 
a = p = 1 . 1 mrad 
T e = 0.2 sec 
F = 6.25 frames/sec 
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Thus, the expression for MRTD became: 



MRTD = 4.00 X 10 -4 f/M(f) (°C) (4-6) 

MTF and MRTD values are found in Table 8. A plot of the 
system CTF, MTF and MRTD values are found in Figures 4.9, 
4.10, and 4.11. 
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V. ANALYSIS OF DATA 



A . GENERAL 

As discussed in Chapter IV, the data presented in 
Figure 4.3 and Table 5 showed that the Thermovision-sensed 
temperatures were substantially lower than the actual 
temperatures measured by the thermisters on the four days of 
data recording. The average Thermovision-thermistor 

temperature differences were -4.5 °C, -1.9 °C, -1.1 °C, and 
-4.7 °C respectively. An explanation was needed as to why 
the Thermovision was underestimating target temperatures. 

B. REVIEW OF THERMOVISION'S THERMAL MEASUREMENT PROCESS 

As detailed in Chapter III, the Thermovision computes 
object temperatures based on two equations: 

*i = T a e o I o + T a ( 1-e o ) T s + (3-~T a )I a (3-2) 



and 



I = A [ Cexp ( B/T ) - l] -1 (3-5) 

The first term on the right hand side of equation 3-2 
represents radiation emitted by the target. The middle term 
represents radiation from the target's ambient atmosphere 
that reflects off the target and propagates to the 
Thermovision. The third term represents atmospheric 
radiation received by the Thermovision. 
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In this thesis I s was assumed to equal I a ,i.e. the 
atmosphere was treated as a continuum extending from the 
target's surroundings to the Thermovision. Thus equation 3-2 
becomes : 

Ii = T a e 0 Io + T a (l-e 0 )I a + (l-T a )I a (5-1) 

Based on this assumption, equations 3-2 and 5-1 had a 
shortcoming. First, the middle term on the right hand side 
of these equations had an implicit value of atmospheric 
emittance equal to unity, i.e. the ambient atmosphere was 
treated as a black body. The third term on the right hand 
side of these equations had l-T a as the emittance of the 
atmosphere existing between the Thermovision and target. 
Since the ambient atmosphere and the atmosphere between the 
Thermovision and target were to be treated as one body, both 
portions would have to have the same emittance. This 
emittance would be 1-T a . The original assumption that the 
ambient atmosphere emitted as a black body is invalid. Thus, 
equation 5-1 becomes: 

x i = T a £ o I o + T a (l-T a ) (l-e 0 )I a + (l-T a )I a (5-2) 

This equation was used to recalculate temperatures in the 
temperature distributions; however, the discrepancy was still 
not resolved. The approach at this point was to fit 
empirically the thermistor temperature data to a modified 
version of equation 5-2. This required setting the first T a 
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term of the middle term on the right hand side of equation 
5-2 to one. Based on this change and a combination of like 
terms, equation 5-2 becomes: 

*i = T a e o I o ) ( 2 -Eq ) ^-a (5-3) 

Theoretically, this modification could not be justified 
but it did result in good agreement between Thermovision- 
computed temperatures and thermister measured temperatures. 
The implication of equation 5-3 is that the atmosphere 
exhibited no attenuation of ambient radiation reflected from 
the target. 

C. REVISED TEMPERATURE DISTRIBUTIONS 

Utilizing equations 5-3 and 3-5, temperature 
distributions of the Point Sur were recalculated with a 
Hewlett-Packard Model 97 programmable calculator. The 
revised distributions are enclosed as Figures 5.1 through 
5.4. 

A number of observations about these distributions could 
be made. The warmest part of the ship was usually the smoke 
stack area which had temperatures clearly exceeding 20 °C in 
most cases. The one exception was the starboard side on 6 
May. Element B7 was the warmest part in this case (19.8 °C) . 

Excluding the smoke stack, the warmest part of the ship 
varied from day to day. These areas are indicated below: 
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DATE 


LOCATION (element) 


SIDE 


TEMPERATURE 


6 May 


B7 


Starboard 


19.8 °C 


9 May 


All 


Port 


o 

o 

o 

CM 


11 May 


B8,B9 


Starboard 


o 

o 

00 

o 

(N 


13 May 


C2 


Starboard 


20.7 °C 



The higher temperature associated with row A (just above 
the water line) were most likely due to additional thermal 
radiation from the ship reflecting off the sea surface. The 
higher temperatures associated with row B corresponded to an 
enclosed portion of the ship. Apparently, heat was 
transferred from these enclosed spaces to the outer surface 
of the ship. The higher temperature at C2 was possibly due 
to solar glint. 

The coolest part of the ship also varied on a daily 
basis. These areas are noted below: 



DATE 


LOCATION (element) 


SIDE 


TEMPERATURE 


6 May 


Cl 3 


Starboard 


15 . 2 


°C 


9 May 


A2 


Starboard 


13.8 


°C 


11 May 


A15 


Port 


14.3 


°C 


1 3 May 


A2 


Port 


14 . 8 


°C 


Generally , 


these cooler portions 


were at 


extreme points 



on the ship's superstructure. Exposure to the open air or 
sea surface at these points probably resulted in greater heat 
dissipation than at other locations of the ship. 

A final note on the distributions is that the temperature 
difference between adjacent distribution elements rarely 



45 



exceeded the estimated deviation in temperature due to 
experimental error (1.5 °C) except around the smoke stack; 
however, the range of temperatures within the entire 
distribution, disregarding the smoke stack, clearly exceeded 
this error estimate in all cases. This range was 4.6 “C, 

6.6 °C, 6.5 °C, and 5.9 °C respectively for each of the four 

days . 

D. PRECISION ESTIMATES FOR THERMAL MEASUREMENTS 

The estimate of experimental error was based on the 
differential method using equations 5-1 and 3-5. The error 
in total thermal units due to errors dT a , de G , and dT was: 

dlj_ = 3I/9T a dT a + ai/9E 0 de Q + 9I/9T dT 

= [e 0 I 0 + ( e o -2 ) *a 3 ^ T a + t^a^o + ^ E o 

+ [(I i 2 B/AT 2 ) exp ( B/T ) ] (5-4) 

Subsituting in estimated values dT a = 0.01, de 0 = 0.01, and 
dT = 0.2 °C and typical values for the I terms, the estimate 
of thermal value error was dl - 0.75 t.u. 

The subsequent error in temperature calculations due to 
thermal value error was: 

dT = aT/ai-L dli = [(T 2 /I i B) / (1 - I/A)] dl j_ (5-5) 

Subsituting typical values for temperature and thermal units, 
the estimated error in temperature measurements was 
dT - 1.0 °C. 
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An additional error in temperature values was introduced 
during the process of making the 6 X 15 element temperature 
distributions from images displayed on the Thermovision's 
computer CRT. Depending on the range of the Point Sur, two 
to four CRT pixels were averaged to establish a temperature 
for each element in the 6 X 15 arrays. Temperatures varied 
by as much as 1 °C between CRT pixels. The estimated error 
attributed to this averaging process was 0.5 °C. Thus the 
total estimated temperature error was: 

Total error = Measurement error + Averaging error 
= 1.0 °C + 0.5 °C = 1.5 °C 

E. REVIEW OF REVISED THERMOVISION VS. THERMISTOR 

TEMPERATURE DIFFERENCES 

The recalculated temperatures were also compared to the 
thermistor temperatures. These results are found in Table 9. 
Based on review of Table 9, Thermovision-computed 
temperatures agreed with thermistor temperatures (within 
experimental error) in 20 out of 26 cases during the four 
days of measurements. Six additional cases could not be 
compared because temperatures exceeded the Thermovision's 
thermal range setting. 

In the six cases in which agreement did not exist, few 
conclusive trends about these anomolies could be established. 
These cases involved four different thermistors over the 
course of all four days. One notable observation was that in 
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five cases the Thermovision's predicted temperatures were 
much lower than the thermistor temperatures. Apparently the 
thermal radiation was partially shielded or severely 
attenuated in these cases. 

A point that must be considered is that each thermistor 
was measuring an area of ship surface on the order of one 
square centimeter while each element in the Thermovision's 
temperature distributions corresponded to an area on the 
order of one square meter. Thus, any localized heat sources 
on the ship may have been obscured by more dominant ambient 
conditions . 

F. MTF AND MRTD CURVES 

One observation of the MTF/MRTD plots (Figures 4.10 and 
4.11) was the cutoff frequency, i.e. the spatial frequency at 
which the Thermovision could no longer resolve a target. 
This frequency was 400 cycles/radian; however, this value has 
little utility because it corresponds to an infinite MRTD. A 
more practical frequency was the critical frequency which was 
determined by extending a line through the approximately 
straight portion of the MTF curve to the intercept of the 
frequency axis. This critical frequency was 360 cycles/ 
radian which corresponded to an MRTD of 2.2 °C. As an 
illustrative example: the critical frequency would equate to 
a lateral distance of 1.4m for a target at 1000 m. Thus, 
the Thermovision could resolve target structures to within 
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1.4 m provided that the target to background temperature was 
at least 2.2- “C. 

The MRTD curve was also used to determine the MRTD of an 
element from the 6 X 15 element temperature distributions. 
Recalling that these elements represented 2.75 m width along 
the ship's length, the horizontal spatial frequency of an 
element is 182 cycles/radian at a target range of 1000 m. 
This frequency corresponds to an MRTD of 0.2 °C. An 
important observation is that this MRTD is much less than the 
uncertainty in temperature due to experimental error 
(1.5 °C) . Thus, this pixel size did not limit the 

experiment. However, pixel size could be a factor if the 
MRTD was comparable to the temperature uncertainty. The 
limiting element size in this case would be 1.6 m based on an 
MRTD of 1.5 °C. Thus, the Point Sur with a length of 41.2 m 
could be divided into a maximum of 25 horizontal elements for 
temperature distributions. 

A clarification at this point should be made with regard 
to vertical resolution. The system MRTD curve developed in 
this thesis was based on horizontal spatial frequencies; 
therefore, the analysis of vertical resolution was not 
possible in this thesis. 

The Thermovision’s computer displayed the system's 
7° X 7° field of view with the array of 128 X 128 pixels. 
This pixel size corresponded to a spatial frequency of 525 
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cycles/radian. Thus, the image video would be limited by the 
Thermovision's components preceding the computer display. 

A final observation of the MTF and MRTD plots was that 
these curves were very typical except for a portion of the 
MTF at the higher spatial frequencies. In this region the 
slope of the curve for a typical MTF would be monotonically 
decreasing. This was not the case for the Thermovision's 
MTF; however, this was most likely due to data scatter, 
rather than the Thermovision having atypical performance 
characteristics . 
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VI. CONCLUSIONS AND RECOMMENDATIONS 

A . SUMMARY 

Thermal images of the Point Sur were recorded at two 
locations on Monterey Bay with the AGA Thermovision 780. The 
Thermovision's temperature evaluation program was calibrated 
under laboratory conditions with a black body source set at 
known temperatures. This program was used to establish 
temperature distributions of the Point Sur images. These 
temperature distributions consisted of 6 X 15 element arrays 
which portrayed the ship's superstructure. 

Temperatures from eight elements on these distributions 
were compared to temperatures measured by thermistors. These 
measurement devices had been installed during the recording 
of the ship's images and were located at eight points on the 
Point Sur that corresponded to the eight elements in the 
temperature distributions. A deduced correction was made to 
an equation used in the Thermovision's temperature evaluation 
program, so that Thermovision-computed temperatures agreed 
with thermistor-measured temperatures. 

System CTF, MTF , and MRTD curves were also produced. The 
MRTD curve was used to determine the resolution limit of the 
Thermovision and to compare this limit with the resolution 
limit of the temperature distributions and the 
Thermovision's computer display. 
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B. CONCLUSIONS 



Initial comparison of target temperatures computed by the 
AGA Thermovision 780 and actual temperatures measured by 
thermisters revealed that the Thermovision underestimated 
actual temperatures on the four days on which data were 
recorded. Two modifications were made to an equation used by 
the Thermovision's computer to calculate target temperatures. 
The first modification was made based upon the assumption 
that the target's ambient atmosphere was essentially 
identical with the atmosphere between the target and 
Thermovision. This assumption had the implication that the 
ambient atmosphere's emissivity would be the same as the 
other portion of the atmosphere. The second modification was 
empirically determined and when used in conjunction with the 
first modification, good agreement existed between the 
Thermovision and thermisters to within 1.5 8 C. In the few 
cases in which agreement clearly did not exist, the 
Thermovision underestimated target temperatures. This was 
attributed to the masking of target radiation rather than any 
measurement error introduced by the Thermovision. 

The AGA Thermovision 780 MTF and MRTD curves produced 
from images of aluminum bar targets resembled the curves of a 
typical system depicted in Figure 2.8 except at the higher 
spatial frequencies of the MTF curve. The irregularities in 
this portion of the curve were attributed to data scatter 
associated with the very small target signals in these 
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frequencies. The Thermovision's cutoff frequency was 

determined to be 400 cycles/radian. The critical frequency 
was 360 cycles/radian at 2.2 °C MRTD . 

Another important spatial frequency was that which 

corresponded to the width of an element in the temperature 
distributions. This frequency was 182 cycles/radian at 
0.2 °C MRTD. 

C . RECOMMENDATIONS 

While the Thermovision's temperature prediction 
capabilities proved satisfactory under various atmospheric 
conditions, further testing should be conducted to verify 
that changes to one of the system's temperature prediction 
equations are justified for other conditions. 

The use of 6 X 15 arrays to depict the Point Sur's 
temperature distribution was viable but could be improved. 
First, the number of elements in these arrays could be 
increased in order to reduce the error introduced when 
several CRT pixels are averaged to determine a temperature 
for one array element. The maximum number of horizontal 
pixels would be limited to 25. The Thermovision's MRTD for 
this pixel size would equal the experimental error in 
temperature (1.5 °C) . Second, an automated means of 

extracting temperature information from the Thermovision's 
computer is needed. The manual process of recording 
temperatures from the CRT screen and calculating average 
temperatures is tedious and time consuming. 
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The determination of the system MTF and MRTD curve could 
be improved by using an amplifier to increase the gain of the 
video entering the oscilloscope. This action would increase 
the accuracy of measuring radiation contrast and reduce data 
scatter at the higher spatial frequencies where target 
signals are quite small and difficult to measure. Also, the 
Thermovision should be evaluated using the MRTD test 
described in reference one and the results compared to a MRTD 
curve determined using the method presented in this thesis. 
Finally, MTF and MRTD curves should be developed for vertical 
spatial frequencies, so that the Thermovision's resolution 
characteristics in the vertical direction can be assessed. 

Any further testing with the Thermovision should be done 
with the Thermal Level and Thermal Range settings adjusted so 
that the entire target distribution of a target can be 
determined. One limitation of this thesis has been the lack 
of temperature information around the Point Sur's smoke 
stack. This limitation was due to the Thermal Level being 
set too low, so that the highest temperatures on the ship's 
superstructure could not be measured. 
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APPENDIX A 



FIGURES 




Figure 2.1 Planck's law for spectral radiant emittance 
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Figure 2.2 Spectral distribution curves related to the sun 
(Shaded areas indicate absorption due to 
atmospheric constituants ) 




Figure 2.3 Contributions to background radiation in sky 
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Figure 2.4 Spectral radiance of cumulus cloud 




Figure 2.5 Reflectance and emissivity of water 
average) vs. angle of incidence 



2 to 15 pm 
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MTF and MRTD curves of thermal imaging system 
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Figure 3.1 Thermovision dual scanner with black and white 
chassis monitor 




Figure 3.2 Schematic of total AGA Thermovision 780 system 
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Figure 3.3 Thermal level and thermal range settings vs. 
black body temperature 
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M: Moss Landing Location 
H: Hopkins Marine Station Location 

Figure 4.1 Locations of Thermovision and approximate 
direction of aim 
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Figure 4.2 Thermovision supplemental calibration curve (°C) 
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Figure 4.3 Point Sur temperature distribution ( *C) 
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Figure 4.5 Point Sur temperature distribution ( *C) 
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Figure 4.6 



Point Sur temperature distribution ( # C) 
13 May 1987 
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Figure 4.7 Locations of thermistors 
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Figure 4.8 Bar target pattern perceived by Thermovision 
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Figure 4.9 Thermovision system CTF curve 
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Figure 4.10 Thermovision system MTF curve 
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Figure 4.11 Thermovision system MRTD curve 



72 











CD 


in 




CD "i 


rsr 












. 




t-H 




• 


i 












CD 


C- 


CD 




D ! 


CD 












H 


t-H 






t-H 


t-H 


































CO 


o 


kt cm 




CD 


CD 


CO 1 


1 ^ 












1 — < 








• 


• 








r- 


CD 


A 




CD 


r- 


^ 1 


in 








t-n 


t-H 






t-H 


t-H 


t-H 


! *-H 












d) 
















O 


m 


in 


CO U 




T}* 


CD 


CD 1 


i TT 










. 


t-H 3 






• 










CD 


t- 


t'— 


-P 




CD 


r- 


00 


CD 






t-H 


t-H 


t-H 


CO 




t-H 


t-H 


t-H 


t-H 












u 
















o 


C*~ 


in 


CM 4) 




c- 


CO 


t-H 


t-H 












Q| 






. 


• 








CD 


C- 


00 


6 




CO 


r- 


CD 


00 






t-H 


t-H 


t-H 


<D 




t-H 


t-H 


t-H 


t-H 












«p 
















CD 


CD 


CD 


t-H 


t-H 




CM 


CM 


C- 


00 




* 






. 


• 


i-H 






• 


- 


. 




* 


o 


CO 


00 


00 


* 




t- 


C- 


CD 


CD 






CM 


t-H 


t-H 


t-H 






t-H 


t-H 


t-H 


t-H 








CD 




O 




C*- 


o 


t-H 


in 








• 




H 


. 


. 




. 


. 






CD 


CD 


c- 




m 


in 


r- 


CO 


CD 






t-H 


t-H 


t-H 




t-H 


t-H 


t-H 


t-H 


t-H 


r- 




CD 


c- 


CD 


in 








t-H 


o 


00 


CD 




. 




. 






CD 


. 


. 


. 


. 


. 


. 


in 


in 


in 


m 


CD 


r- 




in 


in 


CD 


t- 


CD 


CD 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 




t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


CD 


r- 


c- 


CO 


t-H 


CD 






C- 


CO 


o 


CD 


t-H 










. 




00 




. 


. 








in 


m 


in 


CD 




t- 




in 


in 


t- 


t- 


CD 


CD 


t-H 


H 


t-H 


t-H 


t-H 


t-H 




t-H 


t-H 


t-H 


t-H 


t-H 


t-H 


t- 






CM 


CM 


CO 




CM 


CD 


CM 


o 


rr 


00 














r- 






. 


. 






in 


in 


CD 


CO 


C*- 


C*- 




CO 


in 


C^ 


t— 


CD 


00 


t-H 


t-H 


t-H 


t-H 


t-H 


t-H 




t“H 


t-H 


t-H 


t-H 


t-H 


t-H 




CD 


CM 


CM 


fH 


CD 








CD 


CM 


t-H 




. 




. 


. 




CO 






. 


. 


. 




in 


CD 


CD 


C- 


t- 








r- 


CD 


CD 




t-H 


tH 


t-H 


t-H 


t-H 








t-H 


t-H 


t-H 




in 




CO 


i C^ 




t-H 


t-H 


CM 


Tj< 


t-H 




. 


. 


. 


. 


in 


• 


. 


. 


• 


. 




CD 


CD 


C- 


C'- 




CO 


00 


00 


00 






t-H 


t-H 


t-H 


t-H 




t-H 


t-H 


t-H 


t-H 


t-H 




o 


CD 


CD 


CM 








in 


00 


in 




• 


• 


• 


• 










• 


• 






CD 


C^ 


00 








m 


C- 


t- 




t-H 


t-H 


t-H 


t-H 








t-H 


t-H 


t-H 




o 


CD 


t- 


in 








CM 


CD 


CD 




. 


. 


• 


. 


CO 








. 


• 




t- 


CD 


C*- 


00 








in 


r- 


c- 




t-H 


t-H 


t-H 


t-H 








t-H 


t-H 


t-H 




CM 


CD 


in 


in 










t-H 


CO 




. 


. 


. 


. 


CM 
















CD 


t- 


00 










CD 


CD 




t-H 


^H 


t-H 


t-H 










t-H 


t-H 




CD 


CD 














CD 


00 




CD 


CD 






t-H 








in 


in 




t-H 


t-H 














t-H 


t-H 


Gh 


w 


Q 


CJ 


OQ 


< 




Ct i 


w 


Q 


o 


CG 


< 



CM 



CO 



lO 



CD 



CD 



CD 



CM 



CO 



in 



Figure 5.1 Revised temperature distribution ( # C) 

6 May 1987 
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Figure 5.2 Revised temperature distribution 

9 May 1987 
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** temperature > 20. 
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Figure 5.3 Revised temperature distribution ( *C) 

11 May 1987 
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Figure 5.4 Revised temperature distribution ( 

13 May 1987 
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APPENDIX B 



TABLES 
TABLE 1 

RADIOSONDE DATA 



DATE 


6 May 


9 May 


11 May 


13 May 


TIME OF LAUNCH 


1905 


0947 


1009 


1008 


TEMPERATURE (°C) 


15 . 31 


11.61 


12.31 


14.01 


DEW POINT (°C) 


12 . 31 


11.51 


12 .11 


13.01 


RELATIVE HUMIDITY (%) 


82 . 20 


99.61 


99.01 


93.91 


PRESSURE (mbar) 


1008.71 


1015.31 


1009.21 


1006.51 


WIND DIRECTION (degrees) 


238 . 31 


335.81 


228 .11 


331.61 


WIND SPEED (m/sec) 


1 .61 


3.51 


1 .91 


8 .11 


ALTITUDE (m) 


29.01 


10.01 


34.01 


55.01 
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TABLE 2 



TEMPERATURES AT SHIP AND TRANSMITTANCE VALUES 



DATE 


6 May 


6 May 


9 May 


9 May 


TIME 


1855-1858 


1859-1903 


1126-1130 


1138-1143 


AIR TEMP. ( °C) 


14.5 


14.5 


12.0 


12.0 


SEA SURF. TEMP. 


(°C) 15.2 


15 . 2 


13 . 2 


13.2 


SHIP RANGE (m) 


800 


800 


1000 


1000 


TRANSMITTANCE 


0.91 


0.89 


0 . 67 


0.62 



DATE 


11 May 


11 May 


13 May 


13 May 


TIME 


0942-0950 


1002-1010 


1518-1523 


1524-1528 


AIR TEMP. ( °C) 


12.1 


12.1 


14.0 


14.0 


SEA SURF. TEMP. 


(°C) 13.4 


13.4 


15.0 


15.0 


SHIP RANGE (m) 


800 


800 


700 


700 


TRANSMITTANCE 


0 . 90 


0.90 


0.93 


0.93 
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TABLE 3 



THERMOVISION SUPPLEMENTAL CALIBRATION DATA 



BLACKBODY 

TEMPERATURE 

(°C) 



THERMOVISION MEASURED 
TEMPERATURE 
( °C) 



TEMPERATURE 

DIFFERENCE 

(°C) 



6.0 ± 0.1 


5.1 ± 0.1 


- 0.9 ± 0.1 


7.0 


6.1 


- 0.9 


8.0 


7 . 2 


- 0.8 


9.0 


15.8 


+ 7.8* 


10.0 


9.1 


- 0.9 


11.0 


10 . 1 


-0.9 


12.0 


11 . 3 


- 0.7 


13.0 


12.5 


-0.5 


14.0 


12.5 


- 1.5* 


15.0 


13.6 


- 1.4* 


16.0 


15.6 


- 0.4 


17.0 


16 . 5 


- 0.5 


18.0 


18.0 


0.0 


19.0 


19.1 


+ 0.1 


20.0 


20.0 


0.0 


21.0 


21.0 


0.0 


22.0 


22 . 1 


+ 0.1 


23.0 


23.4 


+ 0.4 


24.0 


24.4 


+ 0.4 


25.0 


25.4 


+ 0.4 


26.0 


26.6 


+ 0.6 


27.0 


27 . 6 


+ 0.6 


28.0 


28 . 7 


+ 0.7 


29.0 


29.3 


+ 0.3 


30.0 


31.2 


+ 1.2 



* Data point disregarded. 
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TABLE 4 



THERMISTOR DATA (°C) 



DATE 

TIME 

THERMISTOR 1 
THERMISTOR 2 
THERMISTOR 3 
THERMISTOR 4 
THERMISTOR 5 
THERMISTOR 6 
THERMISTOR 7 
THERMISTOR 8 

DATE 

TIME 

THERMISTOR 1 
THERMISTOR 2 
THERMISTOR 3 
THERMISTOR 4 
THERMISTOR 5 
THERMISTOR 6 
THERMISTOR 7 
THERMISTOR 8 



6 May 
1855-1858 
16 . 48±0 . 10 
17 . 39±0 .13 
21 . 11 ± 0 . 22 
17 . 33±0 . 10 

18.7110 . 27 

17 . 2910 . 12 
24.0510 . 34 

21.7610 .13 

11 May 
0942-0950 

15.4110.15 
15.9610.08 

16.0710.06 

15.7910.11 

15 .2310 .19 

16.8910.13 

35 .6210 .23 
39 . 1510.52 



6 May 
1859-1903 

16 . 3110 .13 

16.9710.18 

20 . 9410 .12 

17 .1710 . 10 
18.6410 . 29 

16 . 9410 .15 
24 . 2910 . 55 

21 .7610 . 20 

11 May 
1002-1010 

15.2310.12 

15.8910.12 

15.9210.06 

15.6810.11 
14.9510.22 
16.9010 . 10 
34.5710 . 31 

38 . 2410 . 19 



9 May 
1126-1130 
16 . 5510 . 30 

17 . 7110 .17 

17.4110.11 

17.2510.23 

18 . 2210 . 27 

17.8110.18 

34 .6310 . 24 
33.2210.80 

1 3 May 
1518-1523 

18 . 3910 .16 

19.3110.27 
18 . 8510 .15 

20 .5810 . 20 

26.9710.07 
20.1510.41 

39.5010.18 
39 .6610 .78 



9 May 
1138-1143 

15 . 5810 .16 

17.0010.17 

17.6710.14 

16 .7210 . 15 
18.6710.22 

17.2210.13 
33.5810 . 38 
28 .1610 .40 

1 3 May 
1524-1528 

17.9610.19 

18 .3310 . 19 
20 .2310 . 29 

20 . 1010 . 12 

27 .4210 .15 

18 . 6110 .16 
41 . 3610.59 
37.3410.51 
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TABLE 5 

TEMPERATURE DIFFERENCE (°C) 
(THERMOVISION TEMPERATURES-THERMISTOR TEMPERATURES) 



DATE 




6 


May 


9 


THERMISTOR 


1 


- 


3.4 


- 


THERMISTOR 


2 


- 


3.4 


- 


THERMISTOR 


3 


- 


6 . 2 


- 


THERMISTOR 


4 


- 


3.1 


- 


THERMISTOR 


5 


- 


5 . 3 


- 


THERMISTOR 


6 


- 


3.2 


- 


THERMISTOR 


7 


- 


8 . 2 




THERMISTOR 


8 




3 . 2 





* Not available. 



May 


11 


May 


1 3 May 


1 . 3 


+ 


1.8 


- 3.3 


1 . 3 


- 


1.4 


- 3.2 


3.9 


- 


1 . 2 


- 3.2 


1.8 


- 


1 . 4 


- 3.9 


2.0 


- 


1 . 4 


-12 . 1 


1.0 


- 


2 . 7 


- 2.7 


* 




* 


* 


* 




* 


* 
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TABLE 6 

RADIATION CONTRAST 

(Target Temperature = 25.0 °C, Room Temperature = 18.8 °C) 



SPATIAL FREQUENCY* 


^max 


H min 


C** 






( cycles/radian ) 


(cm) 


(cm) 


(%) 






4.00 


+ 


0.05 


1.40 ± 0.02 


0.20 ± 0.02 


75.0 


± 


1.9 


13.0 


± 


0 . 2 


1.40 


0 . 22 


72.8 


± 


2.0 


25.0 


± 


0.5 


1.40 


0.26 


68.7 


± 


2.2 


36.0 


± 


0 . 9 


1.40 


0 .28 


66.7 


± 


2.5 


46.0 


± 


1 . 3 


1 . 40 


0 .38 


57 . 3 


± 


2 . 3 


56.0 


± 


1 .8 


1 . 40 


0.44 


52.2 


± 


2.2 


63.0 


± 


2.2 


1.40 


0 . 46 


50.5 


± 


2.3 


71.0 


± 


2.8 


1.40 


0 . 48 


48 . 9 


± 


2 . 4 


83.0 


± 


3 . 6 


1 . 40 


0 . 50 


47 . 4 


± 


2 . 6 


100 


± 


5 


1.40 


0.58 


41.4 


± 


2 . 5 


111 


± 


3 


1 . 40 


0.60 


40.0 


± 


1.5 


125 


± 


4 


1.40 


0.62 


38.6 


± 


1 . 6 


136 


± 


3 


1.40 


0 . 64 


37 . 3 


± 


1.2 


143 


± 


5 


1.40 


0.66 


35.9 


± 


1.6 


167 


± 


6 


1.40 


0 . 70 


33.3 


± 


1.5 


188 


± 


5 


1.40 


0.74 


30.8 


± 


1 . 1 


200 


± 


9 


1.40 


0.84 


25.0 


± 


1 . 3 


215 


± 


7 


1 . 40 


0.86 


23 . 9 


± 


1.0 


222 


± 


5 


1.40 


0 . 88 


22 . 8 


± 


0 . 7 


250 


± 


9 


1.40 


0.94 


19.7 


± 


0.9 


271 


± 


10 


1.40 


0.96 


18 . 6 


± 


0 . 8 


286 


± 


9 


1.40 


0.98 


17.6 


± 


0 . 7 


300 


± 


13 


1.40 


1.10 


12.0 


± 


0.6 


313 


± 


11 


1.40 


1 . 20 


7 . 69 


± 


0.33 


333 


± 


12 


1 . 36 


1.18 


7 . 09 


± 


0 . 31 


350 


± 


15 


1 . 30 


1 . 14 


6.56 


± 


0 .33 


375 


± 


16 


0.65 


0.60 


4.00 


± 


0.23 


400 


± 


17 


0.60 


0.60 


0.00 






* df 




(l/2w)dR + 


( R/2w 2 ) dw 










Where < 


df 


is estimated frequency 


error, w is bar 


width. 


R 


is range 


target 


/ 1 


dw = 0.0002 


m, and dR = 


0.01 m . 









** dC - [2(H ma j. H m j_ r 1 ) / ( H max +H min ) 2 ] dH + ( df/f )C 
Where dC is estimated radiation contrast error. 



82 



TABLE 7 



NORMALIZED RADIATION CONTRAST AND MTF 



SPATIAL FREQUENCY NORMALIZED C MTF* 

(cycles/radian) (%) (%) 



4.00 


± 


0.05 


100 


± 


3 


99 . 7 


± 


3 . 


13.0 


± 


0 . 2 


97.1 


± 


2 . 7 


93.9 


± 


2. 


25.0 


± 


0 . 5 


91 . 6 


± 


2.9 


88.0 


± 


2. 


36.0 


± 


0.9 


88 . 9 


± 


3 . 3 


80 . 3 


± 


3 . 


46.0 


± 


1 . 3 


76 . 4 


± 


3 . 1 


69.2 


± 


3 . 


56.0 


± 


1 . 8 


69.6 


± 


2.9 


64.0 


± 


2 . 


63.0 


± 


2.2 


67.4 


± 


3 . 1 


62.1 


± 


3 . 


71.0 


± 


2 . 8 


65.2 


± 


3 . 2 


58.4 


± 


3 . 


83.0 


± 


3 . 6 


63.2 


± 


3 . 5 


56.5 


± 


3 . 


100 


± 


5 


55.2 


± 


3 . 3 


47 .5 


± 


3 . 


111 


± 


3 


53 . 3 


± 


2.0 


44 . 3 


± 


2 . 


125 


+ 


4 


51.5 


± 


2.1 


41 . 8 


± 


2 . 


136 


± 


3 


49 . 7 


± 


1 . 6 


39.0 


± 


1 . 


143 


± 


5 


47 . 9 


± 


2 . 1 


37 . 6 


± 


2 . 


167 


± 


6 


44.4 


± 


2.0 


34 . 9 


± 


2 . 


188 


+ 


5 


41 . 1 


± 


1.5 


32 . 3 


± 


1 . 


200 


± 


9 


33 . 3 


± 


1 . 7 


26 . 2 


± 


1 . 


215 


± 


7 


31 . 9 


± 


1 . 3 


25 . 1 


± 


1 . 


222 


± 


5 


30 . 4 


± 


0 . 9 


23 . 9 


± 


0 . 


250 


± 


9 


26 . 2 


± 


1.2 


20.6 


± 


1 . 


271 


± 


10 


24.9 


± 


1 . 1 


19.6 


± 


1 . 


286 


± 


9 


23 . 5 


± 


0.9 


18 . 5 


± 


0 . 


300 


± 


13 


16.0 


± 


0.8 


12 . 6 


± 


0 . 


313 


± 


11 


10 . 3 


± 


0.4 


8.09 


± 


0 . 


333 


+ 


12 


9.45 


± 


0.40 


7.42 


± 


0 . 


350 


± 


15 


8 . 75 


± 


0.40 


6.87 


± 


0 . 


375 


± 


16 


5 .33 


± 


0 . 30 


4 .19 


± 


0 . 


400 


+ 


17 


0.00 













* dMTF = dC 

Where dMTF is the estimated MTF error. 



0 

7 

9 

3 

1 

9 

1 

2 

5 

3 

0 

1 

6 

1 

0 

5 

7 

3 

9 

2 

1 

9 

8 

40 

40 

40 

30 
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TABLE 8 



MTF AND MRTD 



SPATIAL FREQUENCY MTF MRTD* 

(cycles/radian) (%) (°C) 



4.00 


± 


0 . 05 


99.7 


± 


3.0 


1 . 60 


± 


0 . 07 


X 


10" 


13.0 


± 


0 . 2 


93.9 


± 


2.7 


5.55 


± 


0 . 24 


X 


10" 


25.0 


± 


0 . 5 


88.0 


± 


2.9 


1 . 14 


± 


0.06 


X 


10" 


36.0 


± 


0 . 9 


80 . 3 


± 


3 . 3 


1 . 79 


± 


0 . 12 


X 


10" 


46.0 


± 


1.3 


69.2 


± 


3.1 


2.66 


± 


0.19 


X 


10" 


56.0 


± 


1 . 8 


64.0 


± 


2.9 


3 . 50 


± 


0.27 


X 


10" 


63.0 


± 


2.2 


62.1 


± 


3.1 


4.08 


± 


0 . 30 


X 


10" 


71.0 


± 


2.8 


58 . 4 


± 


3.2 


4.86 


± 


0.46 


X 


10" 


83.0 


± 


3 . 6 


56.5 


± 


3.5 


5.88 


± 


0.61 


X 


10" 


100 


± 


5 


47 . 5 


± 


3 . 3 


8.42 


± 


1 .01 


X 


10" 


111 


± 


3 


44 . 3 


± 


2.0 


1.00 


± 


0.07 


X 


10" 


125 


± 


4 


41 . 8 


± 


2.1 


1 . 20 


± 


0.10 


X 


10" 


136 


± 


3 


39.0 


± 


1 . 6 


1 . 39 


± 


0.09 


X 


10" 


143 


± 


5 


37.6 


± 


2 . 1 


1 .52 


± 


0 .14 


X 


10" 


167 


± 


6 


34.9 


± 


2.0 


1 . 91 


± 


0 .18 


X 


10" 


188 


± 


5 


32.3 


± 


1.5 


2.33 


± 


0 . 17 


X 


10" 


200 


± 


9 


26.2 


± 


1 . 7 


3.05 


± 


0 . 34 


X 


10" 


215 


± 


7 


25 . 1 


± 


1 . 3 


3.43 


± 


0.29 


X 


10" 


222 


± 


5 


23.9 


± 


0 . 9 


3 .72 


± 


0 . 22 


X 


10" 


250 


± 


9 


20 . 6 


± 


1 . 2 


4 . 85 


± 


0 . 46 


X 


10" 


271 


± 


10 


19 . 6 


± 


1 . 1 


5.53 


± 


0 . 51 


X 


10" 


286 


± 


9 


18 . 5 


± 


0 . 9 


6.18 


± 


0 . 50 


X 


10" 


300 


± 


13 


12.6 


± 


0 . 8 


9 . 52 


± 


1.02 


X 


10" 


313 


± 


11 


8.09 


± 


0.40 


1 . 56 


± 


0.13 






333 


± 


12 


7.42 


± 


0.40 


1.80 


± 


0 .16 






350 


± 


15 


6 . 87 


± 


0.40 


2.04 


± 


0.21 






375 


± 


16 


4 .19 


± 


0 . 30 


3 . 58 


± 


0 . 31 






400 


± 


17 


0.00 

















* dMRTD = (4.00 X 10~ 4 ) [df/MTF + ( f /MTF 2 ) dMTF ] 
where dMRTD is the estimated MRTD error. 



3 

3 

2 

2 

2 

2 

2 

2 

2 

2 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 

1 
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TABLE 9 

REVISED TEMPERATURE DIFFERENCE ( 8 C) 
(THERMOVISION TEMPERATURES - THERMISTOR TEMPERATURES) 



DATE 




6 


May 


9 May 


11 


May 


13 


May 


THERMISTOR 


1 


- 


0 . 4 


- 0.3 


+ 


2.0 


- 


0 . 9 


THERMISTOR 


2 


- 


1.0 


- 0.1 


+ 


0.9 


- 


O 

00 


THERMISTOR 


3 


- 


3.9 


- 2.9 


+ 


0 . 9 


- 


0.8 


THERMISTOR 


4 


- 


0 . 7 


- 0.4 


+ 


0 . 7 


- 


1.4 


THERMISTOR 


5 


- 


2.9 


- 0.8 


+ 


0 . 9 


- 


9.7 


THERMISTOR 


6 


- 


0.8 


- 0.2 


- 


0 . 3 


- 


0 . 3 


THERMISTOR 


7 


- 


6.3 


* 




* 




* 


THERMISTOR 


8 


_ 


0.9 


* 




* 




* 



* Not available. 
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